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Abstract 

We study intermittent maps from the point of view of metastability. Small 
neighbourhoods of an intermittent fixed point and their complements form pairs of 
almost-invariant sets. Treating the small neighbourhood as a hole, we first show 
that the absolutely continuous conditional invariant measures (ACCIMs) converge 
to the ACIM as the length of the small neighbourhood shrinks to zero. We then 
quantify how the escape dynamics from these almost-invariant sets are connected 
with the second eigenfunctions of Perron-Frobenius (transfer) operators when a 
small perturbation is applied near the intermittent fixed point. In particular, we 
describe precisely the scaling of the second eigenvalue with the perturbation size, 
provide upper and lower bounds, and demonstrate convergence of the positive 
part of the second eigenfunction to the ACIM as the perturbation goes to zero. 
This perturbation and associated eigenvalue scalings and convergence results are all 
compatible with Ulam's method and provide a formal explanation for the numerical 
behaviour of Ulam's method in this nonuniformly hyperbolic setting. The main 
results of the paper are illustrated with numerical computations. 

Keywords: Escape rate, Pomeau-Manneville maps, Perron-Frobenius operators, Ulam's 
method, Young towers 

2010 MSG numbers: 37E05, 37M25, 37D25 

1 Introduction 

For maps T : X — )■ X with good expansion properties, quantities such as the rate of decay 
of correlations and the rate of escape into a hole may be read off the spectrum of the 
corresponding Perron-Frobenius (or transfer) operator V : {B, || ■ ||) O acting on a suitable 
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Banach space {B, || ■ || ) of regular observables. The rate of decay of correlations is controlled 
by the gap in the spectrum of the Perron-Frobenius operator between the unit circle {z G 
C : |z| = 1} and the rest of the spectrum of V. Of particular interest for metastability is 
the situation where this gap is between the eigenvalue 1, with corresponding eigenfunction 
(p, the density of an absolutely continuous invariant measure (ACIM), and a second real 
eigenvalue A2 G (0, 1) with corresponding eigenfunction ijj Esj]. 

The notions of metastability and almost-invariant sets for closed systems are fre- 
quently explored via sub-unit eigenvalues of Perron-Frobenius operators, and have been 



used in a var iety of real applications [16|, |2l|, |30|, l3l| ; also, extensions to random dynamical 



systems [28|, l29j have been developed recently. 

If a small hole H is excised from the domain X, as in early work of Pianigiani and 
Yorke jiij, one defines a conditional Perron-Frobenius operator Vx\H{f) '■= Xx\h ■ Vif ■ 
Xx\h)- If T has good expansion properties, a conditional ACIM (ACCIM) may exist |46j, 
with density ipn satisfying Vx\h'^h = ^h'^h, < 1- The density ipn describes the 

distribution of trajectories that have not yet fallen i nto the hole and the geometric rate 
at which trajectories of T fall into H is given by \h 



19, 46 



In recent years there has been renewed interest in open systems and escape rates 



stimulated in part by a series of papers by Collet et al. |l0l-ll3l| and, more recently, by 



a survey paper of Demers and Young [19|]. Much work has been carried out to show 
the existence of physically relevant conditionally invariant measures in various settings. 



First results [13|, |46|] were obtained for maps that are Markov when the hole is removed. 



Shortly thereafter, similar results emerged for non-Markov expanding maps of the interval 
l40l . |49|, unimodal maps |l7l . l33| . CoUet-Eckmann maps [7| and hyperbolic maps on 
^ '20|. For sufficient existence conditions in a general class of dynamical systems, see 



The convergence of conditionally invariant measures to the invariant measure of 
the closed system as the hole closes has been investigated by Demers et al. 0, [13, IS 



using Young towers [50(]. Keller and Liverani [37| provide formulae for the variation of 
the escape rate with the size of small holes. 

Another area of recent study is the connection between metastability and escape rates. 
For example [ill show that under fairly general conditions, the existence of a "second 
eigenvalue" A2 G (0, 1) allows one to find a set A such that both the escape rate into A 
and from A are slower than — log A2. Specialising to uniformly expanding interval maps, 
Keller and Liverani jsTj consider a T with two invariant disjoint intervals perturbed by 
a stochastic kernel; a formula is developed for lime_j.o(l — A2,e)/e as the amplitude of the 
kernel shrinks. In 3J] a similar setup is considered (with deterministic perturbations) 
and it is shown that the first and second eigenf unctions, and ip, of the transfer operator 
for a metastable system can be written asymptotically as a convex combination oiipu and 
^x\H as the metastability becomes complete stability (no communication). Dolgopyat 
and Wright [23] show that in the setting of 3^ with m disjoint invariant intervals each 
supporting a mixing ACIM, one can construct an m-state Markov chain with state-to- 
state jump time distributions that precisely estimate the interval-to-interval jump time 
distributions of small perturbations of the original non-ergodic map. 

So far, the use of Perron-Frobenius operators to study the rate of decay of correlations, 
metastability, the rate of escape, and the relationships between them has been applied 
exclusively to the standard test-bed of uniformly expanding interval maps or those that 
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can be modelled by Young towers with exponentially decaying tail^. In this work we 
move beyond uniform expansivity to study in detail the behaviour of these quantities 
for a nonuniformly expanding map with a finite ACIM, the Pomeau-Manneville (PM) 
map T : [0, 1] — )■ [0, 1] (see equation (T) in Section [2]). The natural Young towers have 
polynomial, rather than exponential, tails and it is well known that the Perron- Frobenius 
operator lacks a spectral gap on the Banach spaces of functions of bounded variation or 
Holder observables typically used in the expanding setting. The reason for the lack of 
exponential mixing is long escape times from small neighbourhoods of the origin, due to 
the presence of an indifferent fixed point at 0. 

We link quantitatively almost-invariance (or metastability) of small neighbourhoods 
[0, Co] with mixing rates of a closed system and the escape rates out of the two open 
systems T|[o,eo] and T|[eo,i]. Furthermore, we relate the second eigenfunction of the closed 
system and leading eigenfunctions (ACCIMs) of the open systems in the small hole limit. 
We connect these ideas with numerical approximations of the ACIMs and ACCIMs via 



Ulam's method 48|. In particular, a small amount of averaging in the neighbourhood 



of can induce a spectral gap, stabilizing the eigenvalue 1 (of V). Although the form of 
averaging we introduce looks artificial the resulting analysis explains a phenomenon of 
spectral gap scaling which is evident when applying Ulam's approximation scheme 48|] to 
the calculation of the ACIM for T. These results (Theorems 15.41 and 15. 6p may be viewed 
as intermittent counterparts to spectral stability results for uniformly expanding maps in 
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We now give a brief indication of some of our main results. 



Two open systems. First, we consider two open systems: the dynamics on a small 
neighbourhood of the origin [0, eo] (with trajectories leaving once they expand outside 
this small interval), and the dynamics on the rest of the domain [eo, 1] (with trajectories 
leaving once they enter a small neighbourhood of the origin). 

1. Dynamics on [0,eo].' There is no ACCIM for the open dynamics on [0,eo]. How- 
ever, a computation of escape rate with respect to a variety of natural probability 
measures with support on [0, eo] (which effectively smooth away the nonexpansion) 
yield^ 1 — \h ^ (see Section [3]). 

2. Dynamics on [eo, I].' The open dynamics on [eo, 1] does possess an ACCIM (Theo- 
rem H^]), which is unique in a certain set of regular measures (Corollary 14.31) where 
the corresponding escape rate is — log ~ 1 — Aj:/ ~ eo (Corollary 14.41) . and is the 
same as for Lebesgue measure. Furthermore, we show that the ACCIM of this open 
dynamics converges in to the ACIM of the closed dynamics on [0, 1] as eo ^ 
(Theorem 14.71) . 



The closed system. Second, in Section |3l we show that gluing together the two open 
systems discussed above to form a simple two-state Markov chain predicts a mixing rate 

^For example, certain logistic maps (l?] ]. 

^We use the following notation: f{x) = 0{g{x)) as x ^ a if and only if there exist positive real 
numbers S and M such that |/(a;)| < M\g{x)\ whenever jcc — a| < 6. We write f{x) ^ g{x) if f{x) = 
0{g{x)) and g{x) = 0{f{x)). 
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or "second eigenvalue" whose gap from 1 scales like the slower of the above two escape 
rates, namely eg. 

We show why this is indeed the case, by explicitly constructing a closed system which 
includes a smoothed expansion near the origin that mimics the effective smoothing ex- 
perienced when computing the escape dynamics on [0,eo]. We prove the existence of a 
second eigenfunction for the transfer operator of this system and upper and lower bounds 
for the second eigenvalue (Theorem 15.41) . 

Our method of proof utilises a novel construction combining two Young towers; one 
with positive mass and the other with negative mass. This construction elucidates the 
connection between the second eigenvalue of this closed system and the two escape rates of 
the two associated open systems, and between the structure of the second eigenfunction 
of the closed system and the ACCIM on [eo,l]. Theorem 15.61 explicitlv describes the 
convergence properties of the ACCIM on [eo, 1] and the second eigenfunction when they 
are mapped down from the tower to the unit interval. 



Application to numerics. A major motivation of the form of smoothing in the neigh- 
bourhood of the origin for the PM map is that our theory is applicable to Ulam's method 
of numerically approximating the Perron-Frobenius operator (sometimes called coarse- 
graining). Ulam's method is the unique Galerkin-projection method that can also be 
interpreted as a small random perturbation of the original dynamics (see Theorem 3.7 
[iijl). In the uniformly hyperbolic setting there is a long hi story of results showing con- 
vergence of Ulam's method for ACIMs jii, 24, 3^, 3^, 42, 43|, and isolated spectral values 
and their eigenfunctions 0, [g], 0, 3^. Applications of Ulam's method to open systems 
for the approximation of escape rates and ACCIMs include Recent work i^] has 

shown that Ulam's method converges for certain nonuniformly expanding maps. 

Ulam's method has been used extensively in a non-rigorous way to detect metastable 
and almost-invariant sets for nonuniformly hyperbolic, multidimensional systems. The 
results obtained are typically very good, when compared with objective measures based 



on physics [30|, IMJ , even though the outer spectral values used can be considered spurious 
and due to discretisation effects. 

We show that our theoretical results explain the behaviour of Ulam's method when 
applied to the PM map, including (i) the estimates of escape rate for the open dynamics 
on [eo,l], and (ii) the scaling of the "spurious" discretisation-induced second largest 
eigenvalue of the Ulam matrix for the closed system with the resolution of the projection 
elements. This present work is a first step in explaining the apparent success of Ulam's 
method on difficult real-world problems. 



Outline of the paper. An outline of the paper is as follows. In Section [2] we describe 
the family of PM maps and known results concerning the ACIMs, and introduce condi- 
tional Perron-Frobenius operators and ACCIMs. In Section [3] we introduce our two-state 
metastable Markov chain model, discuss the two open systems formed by partitioning 
[0,1] into [0, eo] and [eo, 1] and state most of our main results. Section H] contains the 
Young tower setup and proof of the existence of an ACCIM on the larger domain. Sec- 
tion [5] contains the dual positive- negative Young tower construction and proof of existence 
of a second eigenfunction for the closed system with bounds on the location of the second 
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eigenvalue. Finally, Section [6] describes Ulam's method or coarse-graining and presents 
numerical results illustrating the various scalings with hole size. 



2 Invariant measures for Pomeau-Manneville maps 
and conditionally invariant measures 

Let leb denote Lebesgue measure on [0, 1]. Integrals with respect to leb will sometimes 
be written J fd{leb) = J f dx. Let T : [0, 1] O be a Pomeau-Manneville map 4l|, |47l 
which near has the form 

T(x) = x + c«x^+° + c/(x) (T) 

where g is and the derivative g'{x) = o(x") (in conventional little-o notatioiJl). Sup- 
pose also that T has two branches, and let xq be the breakpoint such that T is 1-1 and 
onto (0, 1) on both (0,Xo) and (xq, 1). We suppose also that T is on both (xq, 1) and 
(e, Xo) for every e > and that T' > 1 on both (0,Xo) and (xq, 1). 

When a G (0, 1) these maps support a unique absolutely continuous invariant (prob- 
ability) measure fx (ACIM) (obtainable from a first return ma p l45l| '). the dynamics of 
{T, fj.) is exact, and T exhibits polynomial decay of correlations 50!] with rate 0{k^~^^°'). 
The density ^ = ^ has a singularity at (of power law type with exponent —a), and 
arises as a fixed point of the Perron-Frobenius (transfer) operator V for T. The slow 
decay of correlations occurs because typical orbits of T require anomalously long times 
to escape from the neighbourhood of 0. 



For small eo the set [0, eo) is almost-invariant [15|, |27|]. Often, almost-invariant sets are 
associated with isolated eigenvalues outside the essential spectrum of V [ij] . However, for 
Pomeau-Manneville type map^, the eigenvalue 1 corresponding to the invariant density is 
not isolated from the essential spectrum of V on any reasonable subspace of L^, leaving 
no room for isolated "second" eigenvalues. We will see in Sections [5] and [6] that the 
situation is different for certain approximations to V. 

Escape rates, conditionally invariant measures and conditional 
Perron-Frobenius operators 



We use the following standard concepts (see for example [19|, |46(|). Let (X^) be a 
Borel space and let T : X — )■ X be a measurable, non-singular transformatioiu. For a 
measurable set A let A^'''^ be the /c-step survivor: 

A^''^ ■={xeA : r (x) eA, i = l,...,k}. 

The escape rate from A (or into H := X \ A) , with respect to the measure u is given by 
— log A where log A = limfc_j,po ^°s'^(^^' ^) . in absence of ambiguity we may refer to A itself 
as the geometric escape ratqj, or escape rate. 



3That is lim^r^o ^ = 0. 

■^Indeed any expanding maps with indifferent periodic points. 
^That is, v oT~^ is absolutely continuous with respect to i^. 

is usually called the eigenvalue of A for reasons that will become apparent in the next paragraph. 
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Measures that occur naturally in this setting are those measures, say ua , that satisfy 



i^A o {T\a) = Xl^A- 

Note that ua need be defined only on A but we may extend it to all of X. Then UAiA^''^) = 
and the corresponding escape rate is — log A. These are called conditionally invariant 
measures. 

Recall that the Perron-Frohenius operator for T acts on integrable functions ip G 
L^(X, z/) according to 

(where {ipu){E) = J ipxEdu and ^ is the usual Radon- Nikodym derivative). Non- 
negative, normalised fixed points of V are the densities of absolutely continuous invariant 
probability measures (ACIMs). The conditional Perron- Frobenius operator (with respect 
to a measurable set A) is defined by Va'{' = XAV^ip xa)- In much the same way, non- 
negative normalised eigenfunctions of Va are densities of absolutely continuous condi- 
tionally invariant probability measures (ACCIMs). That is if Va^ = X(p then z/^ := yji/ 
is conditionally invariant. It is often convenient to define the normalised conditional 
Perron- Frobenius operator Va which acts according to 

~ VaV 

Then nonnegative fixed points of Va are densities of ACCIMs, while the denominator in 
dl]) gives the corresponding escape rate. 



3 Coarse graining, met ast ability, and a two-state model 

As indicated above, the eigenvalue 1 is not isolated on any reasonable subspace of 
L^([0, 1], z/). Nevertheless, approximate, typically compact, versions of V derived from 
small random perturbations or numerical schemes possess a spectral gap. One of the 
main goals in this paper is to explain the variation of this gap with the size of pertur- 
bation, or resolution of numerical scheme. To this end, we fix eo G (0,Xo) and partition 
[0, 1] = 1^0,1 U Je(,,2 where /^q^i = [0, eo] and 1^0,2 = (co; !]• For all formal results we assume 
that eo is a preimage of Xq (the hole [0, eo] is Markov). We begin by summarising known 
facts and some main results for T and small perturbations of T on these two sets. 



Mass of and /,,,2. The ACIM has /i([0, eo]) ~ ej"". 



Conditionally invariant measures and rates of escape from leo.jij = 1)2. We 
will show that 

'^Formally speaking, ii X = AU H (where H is a hole) and T : A ^ X, we have Va '■ L^{A, vIa) — > 
L^{X,i'). The definition of transfer operators between different spaces is is easy: if T : (Y, fj.) — > (X,!^) 
is a non-singular transformation (/i o T^^ <C u), then Pt ■ L^{Y, ji) L^{X, v) is the Frobenius-Perron 
operator defined by = ^((V-'m) ° T^^) 
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• there exists an ACCIM fi^^ on /eg, 2, with eigenvalue Xe^ and density ip^o bounded 
away from zero and infinity (shown in Theorem 14.21 but also a direct consequence 
of 46|, Theorem 2]); 

• — logAep, the rate of escape with respect to the ACCIM0 scales as eo (Corollary 

• as eo — )■ 0, the ACCIM converges to the ACIM fi (Theorem 14. 7p in total variation 
norm, provided the passage is through a sequence of Markov holes; and 

• the convergence of the ACCIM on J^g 2 to the ACIM can be modelled numerically 
via Ulam's method (Section |6]). 

The situation for I^^ i is more problematic as the only recurrent dynamics supported on 
is the fixed point at 0, supporting an invariant 5-measure So- There is no "natural" 
ACCImEI. Moreover, leh oT-\h^^^i) = O(min{eo, fc"^/"}) anc0 o = C(A;i-i/"). 

Since the escape from I^^^i is asymptotically polynomial, there is no conditionally invariant 
measure which defines a "natural" rate of geometric escape from Ieo,i- 

"Second eigenvalue" . As pointed out in Section [21 V has no well defined "second 
eigenvalue" . We show that 

• adding a random perturbation on I^^^i which averages according to the push-forward 
of a uniform density on /eQ^2 H T'^I^^ i induces a second eigenvalue A2 with scaling 
1 — A2 ~ eQ (Theorems 15.41 and I5.6p : and 

• coarse-graining the map T by Ulam discretisation with grain size eo produces a 
second eigenvalue A2 with 1 — A2 ~ eg (see Section [H]). 



^And in fact any absolutely continuous measure whose density is bounded away from zero and infinity: 
Let be such a measure. Then there exist positive constants b, B where bip^^ < </> ^ Bip^^ so that 
limfc^oo l/fclog'/"^(4*o,2) ^ liinfe^oo ^/klogB^^gilj^^''^) = logAe^. The reverse inequality is similar. 

natural (cf [13]) ACCIM n' will be one for which {Vi^^^fxho.i ^ ^ fc 00, where Vi,^^ is 
the normalised conditional Perron-Frobenius operator. Fix 770 G [0,eo); define ry^ = /jq.i r\T~^rik-i for 
each fc > (and similarly for e^). Then 







dx tk Cfe 



Note that — > as fc — 00, {efc — C/t+i} is a decreasing sequence (since T is expanding) and there is a 
constant C such that eu — Cfe+i < C e]^" for all k. Let L be such that < 770. Then eu+L < Vk so 

L 



Ck Vk ^]_y2 _ < 1 i _ ,,^^) < 1 i Cefe^+" ^ 



as fc 00. Hence, the measure with density Vf^^ iXieo.i converges weakly to ^o- 

""^^These facts follow because (in the notation of the previous footnote) there is a constant c such that 
ifc-i/" <ek< cfc-i/" and ^ ~ x"". 



7 



A two-state metastable model 



Our first approximate version of tfie Markov operator P is a crude two-state Markov 
cfiain approximation, wliicfi nevertlieless turns out to be an accurate descriptor of the 
important dynamics. For any probabihty measure m one can construct a 2-state Markov 
chain with transition matrix 



p 



1 



K 1 - 



where {Peo,m)ij = .)^"°'^'' - This Markov chain describes the movement between a 

small neighbourhood of and the rest of the interval. The invariant probability vector is 

/_6eO_ "ep \ 

We can view this two-state model in two ways. First, the numbers 1 



what are sometimes called almost-invariance ratios |26l . l27l | and may be interpreted as 
geometric escape rates from the two "open" states. Second, as a coarse-grained "closed" 
system, the eigenvalues of Peo,m are 1 and 1 — a^g — b^^, and so the scaling of the rate of 
mixing as eo — )■ is determined by whichever of a^p, b^o is approaching most slowly. We 
now check how well this two-state model reflects the statistics of T. 

The numbers determined by the choice of measure m and there are several 

natural choices: 

1. m = leb — this choice requires no special knowledge of the dynamics of T and 
corresponds to Ulam's method (discussed below). One has = ~ eg and 

2. m = /i — although no closed formula is available, the general structure of the 
invariant density is well understood, and a^o, b^^ scale as when m = le¥A 

3. m is some combination of conditionally invariant measures on each I^^j — as dis- 
cussed above, while an ACCIM does exist on Ieo,2, no ACCIM exists for /eg^i. 

Let us concentrate on options 1 and 2 above, for which we obtain well defined scalings 
for and b^Q- 



Mass of Jeo^i and Ieo,2- The mass given to the first state by our two-state model is 
^ ~ Cq"" and matches the ACIM scaling of /i([0, eo]). 

Rates of escape from leoj^j = l^S. The rate of escape from the second state is 
— log(l — 6eo) ~ ^eo ~ ^0 and matches the rate of escape from I^g^2 with respect to the 
ACCIM. The rate of escape from the first state is — log(l — a^^) ~ a^g ~ e^. This 
is effectively the rate experienced for the map T if the mass distribution on I^^^i were 
reinitialised to /i or m at each iteration of T. 

f 1+ x~° dx 

"The density |f ~ x"" so a,„ - '"'X" and ^ eo- 
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Second eigenvalue of Peo,m- The second eigenvalue is 1 — — beg. Since aeg + beg ~ e^, 
it matches the scaling of the second eigenvalue of the perturbed operator. 

Thus, this two-state Markov model captures well (i) the relative mass of I^^^i and 
/go^2, (ii) provides escape rates from the two states consistent with true or perturbed 
escape rates for T, and (iii) captures well the mixing rate of a perturbed version of T. 
These properties, expressed very clearly with only two states, will carry across to matrices 
arising from Ulam approximations of V. 

4 Towers, ACIMs and ACCIMs 

We study the open and closed dynamics of T via a Young tower of returns to an interval 
away from the indifferent fixed point at 0. For the ACIM, the construction is standard. 
For the ACCIM, we puncture the tower for the closed dynamics, and look for a fixed 
point of the conditional Perron- Frob cuius operator on the open tower. Because of the 
normalisation in the operator, the fixed point must have growing mass concentration as 
height increases up the tower; some effort is needed to control this growth. 

4.1 Tower set-up and notation 

Recall that T : [0, 1] O satisfies (T) and has two, onto 1-1 branches, with a discontinuity 
at xq. Set Ao = [xq, 1] and let A be constructed as the tower of first returns to Aq. That 
is, let R{x) — min{n > : T"(x) e (xq, 1)} and the tower has height R{x) — 1. Let 
Ao,i = {x e Ao : R{x) = i}. Then A = {{xj) e Aq x Z+ : i< R{x)} and the base of 
the tower Aq is partitioned as {Ao,j}^i. Because {[0,a;o], [xq, 1]} is a Markov partition 
for T, the first return to Aq from the top levels of the tower completely covers Aq. This 
fact is used in our arguments below. The upper levels of the tower are A^ partitioned 
as The tower is equipped with a natural measure u (e.g. Lebesgue on Aq 

lifted by upwards translation). Since R is integrable with respect to z^Iaqj is a finite 
measure. The tower map is F : A O with F{x,i) = {x,i + 1) if £ < R{x) — 1 and 
F{x, R{x)) = {T^{x), 0) with F : (Ao,i x {R{i)}) Aq being injective and onto. Assume 
that F is non-singular and satisfies a variant of the usual regularity condition for x,y on 
the same level of the tower 

JF{x) 

'JF{y) 

where JF = '^'''^f^ ^he Jacobian derivative of F, /3 G (0, 1), c < cxd and s is the usual 
separation time. That is, s{x,y) counts the number of returns to Aq before x and y are 
in different elements of the partition of Aq. 

4.1.1 Truncation and escape from the tower 

Next, for each n we impose a hole in the tower Hn := Ue>i^i>n+iAe/, that is, Hn consists 
of all elements directly above :— Uj>„+iAo,i. Note that the highest level of A \ if„ is 
n — 1 and 

H^^{A\H^)nF-'H^, (2) 



osoF{x,y) 



(JF) 
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that is consists of all the points that fall into the hole in exactly one iteration of 
the map F (see Figure a)). If V : L^{A) O is the usual (Perron-Frobenius) transfer 
operator on the tower then 

and for each hole if„ the conditional operator Vn '■ L^i^) O is given by 



\JF{y)\' 



4.1.2 Distribution of R on the tower 



The distribution of R on Aq is determined by the exponent a. Define the sequence 
{xn} recursively by x„ = (0,xo) fl T~^x„_i. Standard computations give x„ ~ ri"^/" 
and leb{xn+i, Xn) ~ Thus R{x) = i precisely when T{x) G (xj_i,Xj_2) so that 

leb{x : R{x) = i} = 0{i^^^^^") and the tail of the return time distribution is 



J2 leb{R >^}- Yl 



i>n+l i>n+l 

(giving polynomial decay of correlations with rate 0{n}^^^") by 50,]) and 



uiH'j ~ n-i/". (3) 
4.1.3 Satisfaction of (JF) for a tower modelled on T 

Let T be a Pomeau-Manneville map satisfying (T) and let to(x) = min{n > : T"(x) G 
[xq, 1]} denote the first passage/return time to Aq. In order to choose (3 such that (JF) 
is satisfied, note that standard estimates (see for example jsoj) give a constant Cq such 
that 

/T^o ( r) 

<co\r"ix)-T^°{y)\ 



log 



JT-o{y) 

when x,y are in the same 1-1 branch of T'^o. Choosing /3 < 1 large enough that 
\JT{x) \ /3 > 1 for all x E [xi, 1] ensures that \J(T'^°){x)\ > f3~^ for almost every x G (0, 1]. 
Hence, distances between points are expanded by at least (3~^ on every visit to Aq. If 
s{x,y) = n then x,y lie in the same 1-1 branch of (T'^")" so 

|r"«(x) - T"«(y)| < |(T"")"(a;) - (r"°)"(y)| < 

and (JF) follows. 

4.2 Existence and uniqueness of ACCIM 

In this section we prove the existence and uniqueness of an absolutely continuous condi- 
tionally invariant probability measure of F : A with hole Hn- It should be noted that 



Theorems 1 and 2 by Pianigiani and Yorke [46|] may be directly applied in our setting to 
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show these claims. Nonetheless, we duplicate the analogous results on a suitable tower. 
As well as preparing for the sections that follow, obtaining explicit, uniform bounds on 
the density of the ACCIMs (independently of hole size) is necessary for our estimates 
of escape rates, and for showing convergence of conditionally invariant measures to the 
invariant measure as the hole close^. 

For a fixed constant C > let C=k be a set of regular functions in -L^(A, u) defined as 

C* := G Iv^(A, u) : > 0, Lp{x) < (p{y)e^^ ' '''' for a.e. comparable a;,?/| . 

Above, x, ?/ G A are considered comparable if either they are both in the same cell of the 
partition A^ j, or if they are both in Aq. 

Furthermore for every n G Z+ let C„ C C,,, be a family of regular densities on A \ 
that is 

Cn:= eC^ : / LP dv = 1, p\h^ = ^ \ ■ 



Lemma 4.1 For each B > the set Cf := G : ||v5||l°° < B} is compact in 
L^{A, v). In addition for each n there is a B = B{n) such that Cn C Cf so that each Cn 
is also compact. 

Proof. Define a metric djs on A by dp{x^y) = (3^^^'^'^ (for non-comparable x and y set 
s{x,y) = 0). For a given p E and any e > choose 5 = min(/3, C~Mog(l + e/B)). 
Take any x,y E A such that dj3{x,y) < 5. As 5 < /3 we have s{x,y) > 1 so x and y are 
in the same cell of the partition of A. Then 

\v{x) - p{y)\ < ^(x)|l - e^^^'^^'l < 5(e^^''"^' - 1) < i?(e^^ - 1) < e. 

Thus Cf is equicontinuous. Let {(pn} be a sequence in Cf . There is no loss of gener- 
ality in assuming that each \pn\ < B. Since the metric space (A, c?^) is separable, a 
diagonalisation argument similar to the Arzela-Ascoli theorem establishes the existence 
of a pointwise convergent subsequence {pnj}- Clearly the pointwise limit also belongs 
to Cf . Since z^(A) < oo, Lebesgue's dominated convergence theorem implies that {pn } 
converges in L^{A, u). This establishes compactness. 

Now consider p E Cn- For any A^ j, 1 < £ < i < n by the integral mean value theorem 
there exists x* G Ae^i such that p{x*) = J^^ p du/u^AeA), so 

gC/3 r gC/3 



esssup^ < e^^p{x*) = / p du < 



J^{Ae,i 



c 



Similarly on the base of the tower we obtain esssup^^y^ < ?^ If we choose 



1 1 

B = Bin) := e*" max ( — - — -, max 



^Z/(Ao) ' l<i<i<n ulAi^i] 

then Cn C Cf hence C„ is also compact. □ 



-"^^The ACCIMs we construct have uniformly bounded densities on A, but not when projected back to 
the interval [0, 1]. 
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Theorem 4.2 Let C > c/{l — (3). For each n e Z"*" the normalised conditional operator 
Vn admits a fixed point in 

Proof. Note that C/3 + c < C. First, we will show that VnCn ^ Cn from which a 
standard fixed point argument will follow. Let e C„; it suffices to show that Vn^P £ C*. 
Now let {z,i), {w,i) e Ag^i where 1 < £ < i < n so that both {z,i) and {w,i) have only 
one pre-image of F, namely {z,i — 1) and {w,i — 1). Here, separation time is invariant 
under F so s{{z,i), {w,i)) = s{{z,i — 1), {w,i — 1)). Moreover, JF = 1 on these levels, 
since the translation is straight upwards. Hence 



< 



JF{z,i-l) 



1 

((z.o.C"',^)) 



On the base, the story is different as for each (2;, 0), (w, 0) G Aq there are n pre-images on 
the top levels of the tower. For £ = 0, .... n — 1 let {zn^tj G A.^/+i be such that F{z^, i) = 
{z,0) and similarly for {wi,i). Now we have s{{zi,i), {wi,i)) = s{{z,0), {w,0)) + 1 for 
every £ = 0, • • • , n — 1. Then for any (p E Cn 

{Vn^){z,0) = 2^ 



n-l 



\JF{zeJ)\ 



~2^^\JF(wejf 

(C/3+c)/3''((^>o)'('"'0)) 



Since C/3 + c < C we have Vn<p G C„ for all (/? G C„ and therefore P„C„ C C„. 

It is easy to see that C„ is convex as for any G C„ and x,y E A \ Hn, we have 

X^(x) + (1 - A)0(x) < (X^(y) + (1 - A)0(y))e^'^^'^'^\ 
Moreover, the operator Vn is continuous as V is contractive: 

By the integral mean value theorem and the conditions on (/? G , 

Lpdv < e'" —— — / (p dv. 



v{Hl) Jh. " - v{A, \ HI) 
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Hence 



ifdv = j ipdv + I Lfdu + Lfdu 
A Jh^ J^o\Hl Ja\Ao 



< e —n. — r^TTT / (pdv + I if du 



Vn^di'=:a / Vn^dv. 



Since l/HP^v?!! < a, the normalisation map Vn'-P ^ Vn'-P /\\Vn'{>\\L^ is continuous. Com- 
bining the above results with Lemma l^?T] we see that C„ is a compact, convex set, invariant 
under the continuous map Vn- The Schauder Fixed Point Theorem asserts that Vn has 
a fixed point e C^. □ 

We prove uniqueness of in Cn below. 

Corollary 4.3 Let C satisfy the hypothesis of Theorem \4.2\ For each n G Z+ there is a 
unique ipn ^ C„ such thatVnfn = ^nfn where A„ = llPn'/'nllLi ■ In addition ipn is bounded 
above and below by positive constants. 

Proof. Let ipn € C„ be a fixed point of Vn- Then ipn satisfies Vn'^n = ^nfn where 
An = ll'PnV^nlUi- Now if ess inf(y9„ = then the regularity of (pn ensures that (pnlAn = 
a.e. on some A^^j. Take any x G Aq^. Then = A^v9„(F^(a;)) = = 9?„(x), 

hence (Pn\Aoi = 0. All of Aq is comparable so this forces to vanish on Aq and hence 
on all of A. Clearly this is not possible as cpn is a density so necessarily essinfy9„ > 0. 

Now, since cpn is bounded abovj^ and below by positive constants, — log A„ is the 
Lebesgue escape rate into the hole Hn so A„ is unique (c/. footnote [8]) . 



In the proof of uniqueness of we borrow a technique from [46|. Suppose that there 
is another eigenfunction 0„ with the same eigenvalue A„. For any s G M we are able to 
construct another eigenfunction fg of Vn, namely 

fs ■= Sfn + (1 - S)0„. 

Let 0" > 1 be the largest real number so that essinf/^ > for all s G (l,cr]. Then 
necessarily essinf/o- = and f„ = lim^-i^o- fs & Cn- We have already seen in the first part 
of the proof that this cannot be, hence is unique. □ 

Corollary 4.4 Let C be such that Theorem \4.2\ holds and let ipn and A„ be as in Corollary 
\4.3\ For each n G Z+ let fin be the measure on A with density ipn = dfin/du. Then fin is 
an absolutely continuous conditionally invariant probability measure for the open system 
with hole Hn. In particular A„ = 1 — fin{II^). 

Proof. It is a well known result that nonnegative normalised eigenfunctions of the 
Perron-Frobenius operator are densities of absolutely continuous conditionally invariant 
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Note that G C„ and see Lemma 
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probability measures (see e.g. 19|, |46|). As fin is conditionally invariant with eigenvalue 
A„ we have 

Xn = fin{F-\A\Hn)\Hn) 

= fln{F-\A) \ Hn) - fln{F-\Hn) \Hn) = l- flniK). 

□ 



4.3 Convergence of ACCIM to the ACIM of the closed system 

Lemma 4.5 Let G C„ he as in Theorem \4.S\ There exist positive constants a and b 
(independent of n) such that essinfA\//„v^n > and esssup^\^fj^(fn < b for all n G Z+. 

Proof. Fix n G Z+ and let ipn,i = V^nUo,; for 1 < z < n and (pn,n+i = V^nl/fi- First we 
approximate a lower bound of and then obtain a uniform upper bound for A~". We 
begin with the result of Corollary 14.41 to obtain 

A„ = 1 - / ^„d.>l- > 1 - 

where the first inequality above is a consequence of the integral mean value theorem and 
the property that ip{x) < e'^'ip{y) for all x,y E Aq. Now using the fact that i^{H^) ■ n < 
z/(A) we obtain 

n ■ lyyAo) n 

for some constant C := e'"z/(A)/z/(Ao) independent of n. Next choose n* G Z+ so that 
C /n < 1/2 for all n > n*. By the mean value theorem there is a constant C" such that 
log(l — C /n) > —C"/n for all n > n* and hence 



n ^ 

Using the bound on A^" and the fact that is an eigenvector of norm 1 we obtain 

n i 

1= llV^nlUi = ||y^n,n+l||Li + K^^'^^ \\Vn,i\\ 

i=l j=l 

(n i 
Ki/^) + ^5^A;(^-i)KAo.) 
i=i j=i 

< (essinfAo¥^n)e^ l^iK) + ^^"^(^o.i 
\ i=i 

oo 



< (essinfAoV2n.)e^"^'^" ■ ^{A 

= (essinfAo¥^„)e^+'^'V(A). 



04) 

i=l 
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Hence we have a uniform lower bound on on and therefore on A\Hn for all n > n*: 

essM A\H„<fn = essinfAoV^n > (^e'^^'^'V(A) j > 0. 
With dl]) we are also able estimate an upper bound of ipn, for n > n*: 

esssupA\^/„¥5n < A-"esssupAo¥5n < e^"^-r-^ / <^ndu < -YiTT- 

i^y^o) Jao i^y^o) 



Since n* is finite, we conclude that there exist constants a > and b > such that 
essinf(y9„ > a and esssup(y9„ < b for all n > 1. □ 

Corollary 4.6 Let the hypotheses of Theorem \4.^ hold. There are constants a and b 
(independent of n) such that 

log Xn 
a < lim — — < b. 

Proof. Using lim^,^! ^zf = 1 and 1 — A„ = J^i v^n dv in conjunction with the result 
of Lemma 14.51 proves the claim. □ 

Theorem 4.7 For every positive integer n let G C„ and Xn < 1 he as in Corollary 

\4-'J\ Then — )■ ^, where (p is the density of the unique absolutely continuous invariant 
probability measure fi of the closed system F : A O. 

Proof. The result of Lemma [4.51 ensures that all are elements of C^, which, as seen 
in Lemma Kl\ is compact. Hence a subsequence of {'Pn}, say {fm} converges to some 
density ip'. Let {fim} and fi' be the corresponding measures. Then for any measurable 
A C A we have 

/i'(F-M)= lim/i„,(F-iA) 

< hm f^nAF~\A\HnJ) 

= lim Xn,finX^\HnJ 

= lim Xn^finA^) = f^'i^)- 

But fi' o F~^ < fi' is possible only if /i' is invariant, therefore fi' = fi and (p' = ip almost 
everywhere. Hence (pn ^ <p itl L^{A, u) as required. □ 



5 On second eigenfunctions 

The model in the previous section allowed the construction of ACCIMs for certain towers 
with holes. We now apply these ideas to construct a second eigenfunction for a small 
perturbation of the Perron-Frobenius operator for the PM map. The new family of 
towers will be parameterized by e. In connection with the previous section, e will be 
approximately the size of the preimage of the hole Hn- Now, orbits dropping into Hn are 
not lost to the system, but return. The idea is to model the effect of a small random 
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perturbation, analogous to the Ulam discretization; we will use a tower with two bases A'^. 
One piece is essentially A \ Hn, modelling the dynamics of T until leakage into [0, C(e)], 
the other piece will model a smoothed version of the dynamics of T returning from the 
neighbourhood of 0. It should be mentioned that in this section, to distinguish between 
Perron-Frobenius operators of various maps, we may write the corresponding map in 
the subscript (e.g. Vt)- This should not be confused with our previous notation for 
representing the conditional Perron-Frobenius operator. 



5.1 Construction of a revised tower 

Let A be a tower with a polynomial tail as above: 

z/(Ao,) ~ r^-V". (5) 

Next, fix e > and choose n := n{e) to be minimal such that z/{Uj>„+iAo,i} < e. Then 
there is a constant di such that 

n < die-". (6) 
Fix ei = z/(Uj>„+iAo,i) and €2 = z/(Ao,„+i). Then there are constants d2,d3 such that 

d2-<e2<ds-. (7) 
n n 

The first tower is A*"'"*" = A \ Hn where = Uj>„+i U^>i A^^j. Thus, Aq is the base of 
the tower A*"'"*". Put ^Qn^i = Uj>„+iAo,i (the pre-image of the hole Hn) and let {Ao,j}"=i 
(same notation as in Section HTTj) define the remainder of the partition of Aq"*". The base 
of the second tower A^'~ is identified with Uj>„+iAi_j. Let l : (Uj>„_|.iAoi) — )■ [0,ei) 
be a Lebesgue measure preserving bijection such that t(Uj>„+i Ao,i) = [0, ei — £2) and 
t(Ao,n+i) = [ei — ^2, ei)- Put Ag~ = i(Uj>„+iAo,i) x {0}. We will replace the action of F 
on H^ = Ui>„+iAo,i by l; instead of direct upwards translation, H^ is slid out of the way 
and regarded as the base of a separate tower. Partition 

where A^'" = [0, ei - 62) and A'f^~ = [ei - 62, ei). 
The height function is defined as 



h'{x) 



k X G Aq ^ , k < n 

1 X E Aq'„_,_;^ 

1 X G Aq'i 



n X e Aq] 



Then A^ := A^-" U A^'+ where A^-* = {{xj) : x E A^'^i < h'{x)} for * G {+,-}. 
Points in the tower will be said to be "in A*^'*" with the obvious meaning. 
The dynamics on A^'"*" are 



F'{x,£) 



F{x,i) if (x,£)GA^'+\A^ii 
(.(a;),0) (x,0)G A^;:+,. 
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Note that when (x,0) G Aq'+^i, F%x,0) e Aq' . On A^'" put 

F'{x,i) = l (x,i + l) (xj) e Al'-J <n-l 



F"+i(i-i(a;),0) {xJ)eA 



n— 1,71 



Note that when {x,n — 1) G A^'Zi.„, F^{x,n — 1) G Aq'^. 

The picture on A*^'" is different to the usual tower dynamics and has been designed 
to allow the analysis of a small random perturbation of the Pomeau-Manneville map T. 
When the return time to Aq is greater or equal to n (corresponding to arrival by orbits 
of T in a small neighbourhood of 0), the dynamics of T are replaced by weak affine 
expansion near 0. The natural measure z/*^ on A*^ is obtained from Lebesgue measure on 
the interval in the usual way. 

The separation time on A*^ is defined with respect to the partition of Aq^ so the 
satisfaction of (JF) for F*^ is inherited from the corresponding property of F. 



5.2 A modified cone and transfer operator 

Fix constants g,G > and C > c/(l — Let e > and assume the constants di,d2, 
and a are all fixed such that A*" is built with n, ei, €2 satisfying and ([7]). 
(j<:,G,9 consists of functions ip satisfying the conditions (Cl.a)-(C2.b): 

(CI. a) ip > a.e. on A'''+ and J^,,+ ip = \; 

(Cl.b) for almost every z,y G A^'"*", ^p satisfies the regularity condition, 

ip{w) ~ ' 
(Cl.c) letting '^i = {i < n), the growth condition 

||^^||oo<G' (i-f)^'; 

holds. 

(C2.a) < on A^'~ and /^,,_ ^ = ^; 

(C2.b) ip is piecewise constant on the levels of A^'^. Explicitly, 

n-l 
1=1 

and the numbers {a^j^jQ satisfy the growth condition 

\ai>\ < \ao\A~^ 

where 

A = A(e) :=l_2^. 
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(a) 



: : ; : ^: » • 

; ; i ; ; ; Ai,2 ; 

ffi Ao,„ Ao,„-i Ao,2 Ao,i 
h [^o,ln] *H [7„, 1] H 

(b) m m i • 

A"'", : A''+, : : : : : 

9 9 

: : ■ AE,+ : : : 

; • : '-*n-2,n-l 



i \ 1 : \^^\ 

9 9 9 9 9 9 9 9 9 9 

A e,- A e,- A e,+ A£,+ A£.+ Ae,+ A £,+ 

"0,1 ^0,n "0,ri+l "0,7-1 "0,n-l "o,2 "o,l 

H [0,ei] H H [2^0, 7n] <+' hn,!] H 

Figure 1: The two towers: (a) The tower of Section|l]with base identified with [xq, 1]; (b) 
The two towers of Section |5] with bases identified with [0,ei] and [xo,l]. The left hand 
tower contains negative mass and the right hand tower contains positive mass. The map 



18 



Functions in C^'^'^ have zero mean, are negative on A^'~ and positive on A^'+; they 
have good "horizontal" regularity, and a controlled growth rate up the levels of the tower. 
The more flexible growth rate up levels of the tower in A'^'+ is made for convenience. 

Let Vf^ be the Perron- Frobenius operator for F*^ acting on functions in C^'^'^ . Un- 
fortunately, Vp^i^ will typically not be piecewise constant on A^'~. To get around this 
problem, let A}^ act according to 

( 4.- ^c^^- 

A^i^ix) = I uHAln ^ ^ Ao , (8) 

I ip{x) otherwise, 

and put — Vf<^ ■ Notice that is the transfer operator associated with replacing 
F'^[x) with a random perturbation uniformly distributed on Aq~ whenever F'^{x) G Aq~. 
We will show that the normalised version of £^ has a flxed point in C^'*^'^; this is an 
eigenfunction of with eigenvalue bounded below by A(e). We prepare by establishing 
some upper and lower bounds on elements of C^'*^'^. 

Lemma 5.1 Suppose thatO < ip : A^'+ — > R_|_ satisfies ip < B and ip^x) / ip{y) < A 

a.e. if x,y are on the same level of the tower. For all small enough e there is p > 
depending on B and JF hut not e or A such that {C^ip)\^<^^+ > p \\iP\\l^/J^ a.e. 

Proof. There is no loss of generality in assuming HV'IIli = 1. Let R{x) be the first 
return function to [xq, 1] (defined over all of Aq), and leh be Lebesgue measure on the 
interval. Since / R d{leh) < oo there is a finite K such that 'Ylik>K ^^b{x : R{x) > k} < 
2^. Let Tk = UK<e<jAl'j and let i/jx = V^lrx- Note that for each i > K, the £th level of 
the tower 

A^'+ := U^^^+iA;;+ satisfies z/'(A^'+) < leb{x : R{x) > i}; 
in particular, u'^^Tk) < j^- Now estimate 

ipKdu' <B I du' <]-. 

In particular, X\E,+yr^ i' ^ I- Let e be small enough so that n> K. For each < k < K 
let Pk = i/^(A^'++J/zJ(A^'+) and let p = mmo<k<KPk- Then i^^(Ag+J > Pkiy'{Al+) > 
pu^{A^j^'^) so that (using the integral mean value theorem) 

ipdu'>^ [ ipdv\ 



n-1 ^ K-l ^ K-l 



Now, 

""^0 k=0''^k,k+i k=0''^k,k+i k=Q''^'- 

By (JF) and the regularity of V^, for x.y E Aq ^ we have 



£^ij{y) A^^VF^^iy) VF^^iy) 



< Ae' 
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so that 
□ 

Lemma 5.2 Let e and n be as in the definition o/ and suppose that ip G C^'*^'^ (with 
{ae}^~Q as in (C2.h)). There are constants < S, D < oo (independent of e) such that 
\ai\ < D |ao| for each i < n, and |ao| Ci G [5, 1/2]. 

Proof. Since A = 1 — similar to the proof of Lemma I4.5[ A^" < D for some 
D independent of n. Hence, \ai\ < A~^|ao| < -Dlo-ol- Next, the conditions on ip, and 



equation ([7]), show that 

f f d \ 

/ = \ao\ ei + ^ Wkl £2 < |ao| ei 1 + ^D— < \ao\ ti (1 + Dd^). 



^ ■' " k=l \ k=l 

On the other hand, 

^ ^ 1^1 > / m = \ao\ei. 



1 _ 

2 " 

□ 

Lemma 5.3 Let e > be fixed. Each set C''-^'^ is relatively compact in L^{A'', u^). 

Proof. Each function ip & can be decomposed as V^U^- + V'|ae>+- In view of 
(Cl.a-b), for each ip, 2'?/'|^e,+ G C„. Hence, by Lemma Wl\ if {ipk} is a sequence in C'^'^ 
then {ipk\A^^+} has a Cauchy subsequence. By Lemma [521 and (C2.b), if G C'*^'^ and 
z G A*^'" then 

\ip{z)\ < max la^l < D |ao| < 7- — . 

t<n 2ei 

Using this bound, the fact that each V'fcU^.- is piecewise constant (condition (C2.b)) and 
the Heine-Borel theorem, there is a further subsequence {V'fc;} such that {ipki\A^-} is also 
Cauchy. □ 



5.3 Existence of second eigenfunctions and eigenvalue scaling 

Define the normalised operator by 



Theorem 5.4 There are constants g and G (independent of e) such that, C^C^'^'^ C 
(^e,G,3 guffici^nfiy small e > 0. In particular, the operator has an eigenvector 

g ^e,G,g ^j^ff^ associatcd eigenvalue A*^ such that 1 — A*^ G ( — , 



'12 2e2 
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Comments on the proof: Most of the work is in showing that preserves C^'*^'^. 
Mass of ip is exchanged between A'^'"'' and A*^'" in two ways by C^: some fraction of 
"positively signed" mass of ip is transported from A^'"*" to A^'~ — this corresponds to 
points falhng into the hole H^. The total transferred in this way is r»e,+ ih, and it 

^0,n + l 

arrives under one iteration of F on Aq^_,_]^ to Aq'~ (the base of the other part of the 
tower) . Correspondingly, "negatively signed" mass from level A^Ci „ is mapped onto 
Aq^. The proof involves bounding the amount of mass transferred in these ways, and 
controlling the distorting effects on regularity. 

Proof. Let ip G C^'*^'^. The choice of G is made in (G) below, and g is fixed near 
the end of the proof. The conditions (Cl.b) and (C2.b) are projective, so will hold for 
C^ip if they hold for C^ip. Since F'^ acts by upwards translation on A^'~ and acts by 
averaging on Aq~, it is clear that {C^ip)]^^'- is piecewise constant. 

Preservation of the growth condition on A^'^.- We first consider {C'^'ip)\A<^'- (note that 
< 0). Let {ai} be as in (C2.b) and put a'^ = {C'^iI')\a'~- Since the action of F'^ on 
A^'~ is direct upwards translation, and the averaging occurs only over Aq~, for i > 



^1 



\ae^i\ < A ^(A|ao|). 



(9) 



We now show that A oq < 


\a'o\. Since {F^y^A'^^' = 


^0,1 ^0,n+l' 


[ CiP du' 


= [ VF^^dv'= [ 


^du'+ [ 


JaI'- 







tpdu' 



(10) 



in 



{II) 



Now 



ao dx = Xntto = (ei - £2) Oq. 



As for (//), for x,y E Aq'^, 



< e 



so ipi^) ^ 6 



, /^.,+ ^ du' ^ e^' /^,,+ ^ du' _ e^' 1/2 



z/(Ao 



=: G (G) 



where G is independent of x (and e). Since u^{AQ^_^_-^) = ei, (//) < G ei. To complete 
the estimate on (//), use 6 from Lemma \^7I\ and equations ([7]) and l^: 

// < Gei = - — - £2 < Y -7^ ao £2 < y ^ «o £2 < ao £2 

62 \ti J d2 d2 

for sufficiently small e. Hence, by ffTOl) . 



Now 



/ C^dv' = (J) + {II) < ao(ei -£2) + |ao|e2 = ao(ei -2e2) < -|ao|eiA(e). 



lAo'- 



> |ao|A(e). 



(11) 
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In view of ([9]) and ffTTl) . < and the growth condition (C2.b) is satisfied by 

Regularity control on A^'+.- Note that ^XIa^.+ is the identity, so {C''ip)\j\c,+ = ('Pf^V')|a^.+ - 



Now, if w G A^' where £ > then 



since s^{{z,i), {w,i)) = s^{{z,i — 1), {w^i— 1)). The regularity estimate on Aq'*' is more 
comphcated. Write ijj = (f)+ + (p- where 0* = for * G {+, — }. Let z,w G Aq"*" and 

let Zk G A^'^^ ^ (fc = 1, . . . , n) be such that F'^^Zk) = z and let zq G A,^'Zi „ be such that 
F''{zq) = z; use similar notation for {wk}^^Q. Now, 

n 

Vf^^{z) = J2Mzk)/\JF'{zk)\+<P-{zo)/\JF%zo)\ 

k=l 
n 

k=l 
n 

n 

^0+(«;fc)/|JF^M|+0_K)/|Ji^^K)| 
fc=i / 

< e(^''^+'=)^^^'"™' (P^H + (|0_(«;o)|/|Ji^1«^o)|)(C^'/3 + 2c)r(^'-)) 

using the regularity conditions on 0+ = and (JF), the fact that (f>^{zQ) = 0_(wo) 

a„_i < 0, and s^{zk,Wk) = s^F^Zk), F^{wk)) + 1 = s\z,w) + 1. Thus, 



< 



< 



fe=i 



VF^i'iw) ~ 



\ 



VF^i^yw) 



\ 



(HI) 



(12) 



Since C > j^, the constant b := C - (C /3 + c) > 0. 
Claim: When e is small enough, (///) < b. 
Proof of claim: By Lemma 15.21 



\(p-{yo)\ = < D \ao\ < 



D 



while 



> e ' = e 



2ei 
/e6[xo, 1] 



n 



> ^4 — 

£2 ei 
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for a suitable constant d^. Hence < ^- On the other hand, = V'|a^'+ 



is uniformly bounded^^ so by Lemma 15.11 {A = e*"') there is a constant such that 
(Pf^0+)|a^.+ > 4- Hence P^.V^(w) = = Pir.0+(w)-|0_(wo)|/|^^'(wo)| > 



^5 — 77^— and 

_D_ 

(J/J) < (C /3 + 2 c) < 6 

for small enough e. □ 
Now using the claim and ffT^ . 

Vp^i^iw) — ^ ) — 

that is, (Cl.b) is satisfied. 

Control of \\C'^iI)\\li: From the regularity estimates above, when e is sufficiently small 

I Af>- < and {C''ip)\/^t,+ > 0. Since £^ is a Markov operator, it preserves integrals, 
so J C^ip du'^ = and 



||(/:^^)|a^.-IUi = - / C'i^du' = / Ci/jdu' = ||(£^^)|a^,+ ||li 

and all terms in the above expression are equal to Now, 

IIV^U^ - IIli ~ ||('C^V')|a^-||li = total mass exchanged between A^'~ and A"'" 



1^1 + 



< 




< 


ao J 




D 


< 


— e 




2ei 




1 


< 
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(in view of (E]) and ([7])). Putting g := D d^ + 2 diG we have 

= 2 ||(£^^)|a.,-||li = 1 - 2 (||^|a^,-|Ui - ||(/:^V')Ia-IIli) > 1 - -. 

n 

Control of growth up A'^'+.- The fact that {C'^il))\^e,+ < G follows from (G), since = 
1. Higher up the tower, 

= (^f^^)Ia-+ < (1 - g/n)-\VF.nA-^ < (1 - al^V max^, 

since Pi^^t acts on A'^'+ by upwards translation. Since max^E,+ < (1 — g/n)^^^^^'^ G, the 
growth condition (Cl.c) is satisfied. This completes the proof that C^C^'^'^ c C^'^'^ . 



""^^By (Cl.c), '!/'Ia=.+ ^ :— G(l — is bounded uniformly in e by an argument similar to 

the proof of Lemma H3I 
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Finally, the Schauder Fixed Point Theorem in conjunction with Lemma [5.31 gives the 
existence of a fixed point ip"^ of £^ in C^'*^'^. Evidently, C'^ip'^ = \\C'^ip'^\\Li 'ip'^, so tp"^ is 
an eigenvector for and A*^ = The upper bound on A*^ is obtained via the 

estimate (/) and f lTU]) : 

X' ao = {C'r)U- = ° = - m + W) > y = (1 - e^/ei) ao 

(since (J/) > 0); divide by ao < 0. An immediate lower bound on A^ is given by (l — g/n), 
however a more careful estimate using 

(A(e))-i |ao| > \a,\ = = {XT'\{^'r)Lr\ = (A^)"'l^lAr I = (^^)~>oi 

l,n l,n 0,n 

gives the lower bound A*^ > A(e). □ 

Simplifications to the proof of Theorem 15.41 can be made to establish the existence of 
a fixed point of with the properties 

• ip" > a.e. and (p" dv'' = 1 

• V'^Ia^- = X]"=o c^XaJ " ('^^ constants depending on e) 

• '^^Jw) — s*"'^^ for almost every z,w & A^'^ (where C > c/(l — /?)) 

Since = V^S the q are equal to a single constant ancl^ ~ ^. Hence 

J^e- '^^ dv'' = c" (ei + (n — 1) €2) ~ e^~". If v5''|Af.+ is embedded as a function on A^ then 
Lemma 14.11 can be applied (as in Theorem 14.71) to establish convergence to the density 
of the unique ACIM on A. Explicitly, let ip G L^(A, z/) be the density of the (unique) 
ACIM for F. Then 

V^\a^-+ ^ and both \\'P^\a---\\l^a^),\\¥^\a\a--+\\lHA) ^ ^ (13) 

as e — 0. 



5.4 Realisation of the second eigenfunction for the PM map 

When A arises as the first return tower to [xq, 1] for a PM map T, the preimage of the 
hole Hn is = [a;o,7n) (for some 7„ such that T^'jn) = a;„_i — the (n — l)st preimage 
of Xq). There is also 7„+i such that T(7„+i) = x„ and 7„ — 7„+i = 62 while jn — xq = ei. 
The map l translates [a;o,7n) back to [0,ei) with t{x) = x — xq (instead of applying T) 
and the dynamics of F*^ essentially replaces T with a weakly expanding map on a small 
interval [0,eo). To make this more precise, let T[xo,7n) = [0,a;„_i) and put eg = Xn-i 
(note that eo ^ JT{xo) ei ~ e). Now define : A' [0, 1] by 

ihe( «^ _ / ^^(a;) (a;,£) G A^'+ 

^^'^^ - 1 t'+\l-\x)) {x,i) e A^'- ■ 
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By balancing the mass exchange between A^' and A''+. 
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Note that $^(A^;„) = x„_^_i), $'(A^;+) = (for i > 0) and <l>^ : A^' ^ 

[0, eo) is 1-1 and given by T(a; + a;o). By defining 

T-<,) - { ° : 1 1°::"'- <"> 

is continuous (in particular T^(x„) = x„_i = T(x„)), near 0, is a small perturbation 
of T and satisfies (T^)'(O) = S'o := -^^z^- Moreover, T" o = o F*^, so arises as a 
factor of F"^. 

Note: Define $* : A — )■ [0,1] by linie^o '^'^(3;)- This limit makes sense because if 
X G Ae^k then x G A'^'+ whenever e is small enough that n(e) > /c. For all such e, 
<l>^(x) = T^(x). In particular, $*|ae.+ = ^'^|a=.+ for all e and $* o F = T o $*. 

The eigenfunction from Theorem 15.41 does not push down to an eigenfunction of 
the Perron- Frobenius operator for T (or T*^); a small random perturbation is needed. Put 

_ / {fc+i if T{zk) G [0, eo) . . 

\ T(2fc) otherwise, ^ ^ 

where the are i.i.d. random variables on [0, eo) with density function p^i^z) = — J_i 

where T~^^^ : [xo,a;o + ei) [0,eo). Notice that is the push-forward of the uniform 
density on T^ght[0, eo), and is close to constant when e is small (since T is on the right- 
hand branch). Let Vt be the Perron- Frobenius operator for T and let £^ be the (Markov) 
transfer operator associated with the process in f lT^jl l^. Let LI'' := : L^(A^,u'') — )■ 
L^([0, l],/e6) be the Perron- Frobenius operator for and introduce an alternate nota- 
tion for C: C'p, = C'{= A'^ o Vf^). 

Lemma 5.5 In the notation established above, LI^ o Cp^ = o LI''. 

Proof. Let Vp^ be the Perron- Frobenius operator for the system (F", A^), let be 
given by and define A" : L^[0, 1] O by 

j[^f(z) = S P'^^^ fdileh) z G [0,eo) 
1^/(2;) otherwise. 

Then £^ = o Vt- From ([H]), T[x) G [0, eo) if and only if T\x) G [0,eo), so in fact 
li^^ = A^oVT = A^oVT^. 
Next, observe that 

ei ei 

Now suppose that ■?/' is supported on Aq^. Then, by ([8]), 

WoA'^i) = Wx^-^,-— = PI 'ipdu' = p' Wtpdileb) = A' oWip. 



-•^^The small amount of averaging over [0, eg) means that interpolates between Vt, its Ulam ap- 
proximations and the two-state metastable model. 
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If if) is supported on \ Aq then A'^^ip = so 11^ o = TV- Moreover, since 

$'(A' \ Aq = (eo, I), = W ip. From these two cases 

and hence 

o c'p._ = Wo A% o Vf^ = o o Vf^ = A'o Vt^ o = £^ o 

(the third equahty is because $^ o = T*^ o <|>'^). □ 
Theorem 5.6 

1. For any f E L\[0, 1], leh), ^ Vrf as e ^ 0, 

2. C^rp has an eigenvector /" satisfying Crpf^ = and f ^ f* , where f* is the density 
of the unique ACIM for T . 



3. Op has an eigenvector satisfying C'^h'^ = X'^ (where X'^ is as in Theorem \5.4^ , 
and [h']+ ^ \f* as e ^ 0. 

Proof. 

1. First, from the proof of Lemma [5.5[ = A'^ o Vt so for any / G L}, 

\\C-y-VTf\W= r \Cy -VTf\d{leh)<2 H \VTf\dileb) 
Jo Jo 

as eo ^ since Prf E L^{[0, l],leb). 

2. Next, let ip" G L^(A%z/'') be the (unique, normahsed) fixed point of Cp^ and let 
f = nV'- Let v?* G L\A, v) be the (unique) density of the ACIM for (F, A) and let $* 
be the canonical semi-conjugacy from A onto [0, 1]. Then /* := V^*<f* is the density of 
the unique ACIM for T, and by Lemma ESI := = Ti-'Op.^' = £^nV' = Opf. 
Since $*|a^.+ = ^^\a^'+, 



Ili[o,i] 



lir - r IUmo.i] = \\V^4v*\a\a^.+ ) - Tl'i^'W-) + Tl'iiv* - V')\a.+ ] 

< ||p$*(v'*|a\a=.+)IUmo,i] + l|n'(v9'|A^.-)||LMo,i] + l|n'(v3*|A^.+ - v5'|a^.+)||l 
= II¥'*Ia\a'^.+ IIli(a) + IIv^Ia^.- IIli(A') + IIv^*Ia'^.+ - v^1a<='+IIli(a)- 

All these terms approach as e — )■ by (|T3|) . 

3. Finally, we turn to the second eigenfunction. Let g,G be such that Theorem 15.41 
holds and let ip'^ G C^'*^'^ be the fixed point of O. Put h"^ := U'^ip'^. Then, by Lemma [5.51 
and Theorem 15.41 

C'ph' = CfpU'ilj' = U'Op,^' = X' U'l/j' = X' h'. 
By arguments similar to Theorem 14.71 V'^|a'.+ ^-^^ ^ip* as e — )■ 0. Hence 

:= n^(#|A£.+ ) ^f* as e ^ 0. (16) 
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We need to show that [h'^]~^ has this same hmit. This is done by showing that — 
^^||li[o,i] — > as e — 7- 0. To this end, notice that since iP'^\a^- < 0, we have h"^ < and 
in particular [h'^\^ := max{h'^,0} < k'^. 
Claim: lim^^o /o^(^£ ^ [h^'^) d{leh) = 0. 

Remainder of proof of 3., given claim: Using (fTB]l and the claim, let e — )■ in 

\\f*/2-[hT\\L^ < \\f*/2-k'\\L^ + \\k'-[hT\\L^ = \\f*/2-k'\\Li+ !\k^-[h^]^)d{leh). 

Jo 

Proof of claim: Fix rj = e"" and choose J such that rj G (Notice that J = 

C(?7"°).) Since ip" e C''^'^', (Cl.c) gives V'|a'=.+ < G {1 - g/ny < G' for some constant 
G' independent of e. Now [0, r/] fl supp(fc') = (<l>')-i[0, r^) nA'^'+ C ^]=}-i{'^^)~'^[xj,Xj_i) 
so 

/ k' d{leb) < J2 ip' du' < G' = C( J^"^/") = 

•^0 j=j-i A*^)-M^'.,^.-i) e=i j=J 

by On the other hand, 

-/ {h'-k')d{leb) = / -#|A^.-rfi^'< V |a^|z/^(A^';) < 0(Je°) = C(r7~"e"). 



Hence, 



e=n-j 



{h' - k') d{leb) + / k'd{leb)- / A;^(i(/e6) 

> + (17) 

Then, since h'd{leb) = J^.ip'du' = 0, J^\h'\d{leb) > {-h') d{leb) = d{leb) 
has the same lower bound ( 1T7|) . Since t] = e". 



^ |/i^|rf(/e6) > f h'd{leb) 

ri Jr] 



[ ¥d{leb)> [ [h']+ dileb) = - f c/(/e6) > ^-^(e"^^-")) = f k'd{leb)-0{e 
Jo Jo 2 Jo 2 Jq 



a (1— a)^ 



□ 



6 Numerics 

Ulam's method [48,] is a well-known (and effective) method for studying T numerically 
via its Perron- Frobenius operator. For each N E N, partition [0, 1] into sub intervals of 
length and let Bn be the (finite) cr-algebra obtained by taking unions of elements of the 
partition; Sjy-measurable functions are piecewise constant. Let E^r denote the conditional 
expectation operator K[-\B]sf) acting on function in L}{leb). Ulam's method consists in 
replacing V with the finite-rank operator Vn = Eiv o P. The leading eigenvalue of Vn is 
1, and the corresponding eigenvector (fixed point of Vn) is an approximation to a fixed 
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point of V. Surprisingly (given the absence of a spectral gap for P), these fixed points 



converge to ^ as — )■ oo [5l. 144 



Because Eat is a finite rank projection, each Vn can be represented by an x 
matrix P/y. Each P/y is extremely sparse (having 0{N) non-zero entries), and their 
eigenvalues can be found quickly by iterative methods. Because the dynamics of T are 
transitive, each P/v is ergodic, so the eigenvalue 1 has strictly larger modulus than the 
other eigenvalues. Interestingly, the spectral gap (separation of the largest modulus eigen- 
value 1 from the next largest modulus eigenvalue) scales as A^~". 



6.1 Eigenvalue scaling 

The two-state model of Section [3] showed that when the geometric escape rate from the 
set [0, eo] approached zero more slowly than the escape rate from the set [eo, 1], the gap 
from 1 of the second eigenvalue of the two-state Markov chain scaled like the slower escape 
rate from [0,eo]; namely eg. 

We now replace the two-state model with the 'W-state model" arising from Ulam's 
method. The matrix P/v is row-stochastic, representing the transitions of a finite state 
Markov chain whose ith state is identified with the subinterval Jj := [{i — 1)/N,i/N). 
The indifferent fixed point at can be associated naturally with the subinterval Ji = 
[0, 1/A^) ~ [0,eo). The conditional transition probabilities out of state 1 are (PAr)ii = 
leb{Ji n T-^Ji)/leb{Ji) = 1 -NT-\l/N) and (P/v)i2 = 1 - (i'jv)ii- Thus the geometric 
rate of escape from Ji is — log(P/v)ii ^ 1 — (P/v)ii ~ A^~" ~ e^, and this is of the same 
order as previously computed for the two-state model. 

We find numerically that despite increasing the number of states from two to A^, 
the second eigenvalue of our A^-state Ulam matrix retains the scaling predicted by the 
two-state model when eo = 1/A^, namely 1 — A2(A^) ~ A^~"; see Figure [21 




logN 

(a) 




(b) 



Figure 2: (a) Variation of the second eigenvalue of P/v with A^ and a. (b) Slope of line 
of best fit for each a. Note: Computed by eigs in Matlab, with Ulam matrices for the 
PM map 0, Example 3]. 
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Connection with the Perron-Frobenius operators Cfpe and Theorems 15.41 

and 15.61 prove the existence of a second eigenvaluJ^ A'^. The matrix P/v successfully 
reproduces the dynamics responsible for this eigenvalue, and we now explicitly describe 
the connection. Set e = l/(A^T'(a;o )) and choose n, ei,e2 as in ([S])-([7]). Then T[a;o,a;o + 
ei) =: [0,eo) ~ [0, = J\. The averaging created by replacing T by (see ffT^ ) 
or the small random perturbation (ITSl) "almost linearises" T over [0,1/A^], so that P/v 
is a very near approximatioiJ^ to Theorem 15.41 predicts an eigenvalue of A*^" G 
(1 — 2^,1 — ^). Numerical computations with Pjq for a range of N produce second 
eigenvalues within this range. In fact, the upper limit is a very good estimate; see 
Table [H 





1 - A2(Ar) 


e2/ei 


100 


0.069494728128226 


0.060750416292176 


200 


0.047118990434159 


0.042626262679704 


500 


0.028582682402957 


0.026696029895732 


1000 


0.019751285772241 


0.018706181316717 


2000 


0.013727390048589 


0.013165183357731 


5000 


0.008542396305559 


0.008301674655368 


10000 


0.005988977377968 


0.005866565930472 


20000 


0.004208535921532 


0.004150111773511 


50000 


0.002646628586393 


0.002621525600809 



Table 1: Comparison of 1 — A2(A^) computed numerically as a second eigenvalue of the 
NxN Ulam matrix and the corresponding lower bound €2/61 obtained from Theorems 15.41 
and ESI (a = 0.5). 



6.2 Ulam's method and the escape rate from [l/N, 1] 

We conclude with some simple remarks about how to observe the ACCIMs fin, and their 
geometric escape rates, numerically. The open tower systems constructed in Section H] 
exclude a set corresponding to those orbits which land in [0,a;„) under T. Letting 
$* : A [0, 1] be defined by ^*{x,i) = T^{x), the measure /i o (<l>*)~^ is an ACIM for 
T, and /i„ o ($*)~^ is an ACCIM for T\[x„,i] with geometric escape rate 1 — iJin{H\) (see 
Corollary 14. 4p . In fact, fin{Hl) ~ e, where e ~ (x„ is the nth left preimage of the 
discontinuity point xq for T). Hence, if x„ ~ j^, then one expects the escape rate from 
[xn, 1] to scale like 1/A^. Note that with choosing e ~ 1/A^, the corresponding n in the 
constructions of Sections H] and is n ~ e~" ~ A^". Now partition the NxN Ulam 
matrix P/v as 



■ {Pn)ii 


a^ 


b 


po 



^^More precisely, we show that there is another real eigenvalue very close to 1; based on numerical 
computations we conjecture that the eigenvalue A'^ is indeed the second largest real eigenvalue. 

Since T is locally and eo « -^j the transition probabilities {Pn)ii and (P/v) 12 are close to the 
corresponding entries of the matrix for 11'^ o £^ (note that 11'^ — like the perturbation defined by (|15p — 
averages over a small interval [0, ©(e)). Outside [0, eo]j T is uniformly expanding (albeit rather weakly) 
and C° = Vt- 
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where a, b are (A^ — l)-vectors and is an (A^ — 1) x {N — 1) matrix. In fact, is the 
Ulam approximation to the conditional Perron- Frobenius operator X[i/Ar,i]^( ' X[i/Af,i])- 
Figure |3] we present numerical evidence that the leading eigenvalue A°(A^) for P^ has the 
scaling 1 — A°(A^) ~ independently of a. 




Figure 3: Variation of escape rate from [1/A^, 1] with and a. 



Finally, Theorem 14.71 predicts the convergence of the ACCIMs /i„ to the ACCIM 
as the size of the hole Hn shrinks to 0. We illustrate this convergence numerically as 
follows. For a large A^=k (we have used A^* = 10^), form Pjv^ and calculate the leading 



eigenvector. This is a good approximation to the density of the ACIM for T [4^, and 



we use it as a reference measure. Next, for a sequence of smaller A^^ (we used the values 
from the first column of Table [T]) , calculate the leading eigenvector of P^ . Comparing 
the probability measure induced by these eigenvectors with the reference measure from 
the Ulam aproximation Pn, we see good convergence in Figure HI 
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log N 



Figure 4: Total variation norm error between Ulam approximate ACCIMs with a hole 
[0, and the Ulam approximate ACIM with a bin size of 10~^. A range of different 
PM maps are used, {a is order of tangency of the indifferent fixed point). 
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